d-Wave Pairing Correlation in the Two-Dimensional t — J Model 
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The pair-pair correlation function of the two-dimensional t — J model is studied by using the 
power-Lanczos method and an assumption of monotonic behavior. In comparison with the results 
of the ideal Fermi gas, we conclude that the 2D t — J model does not have long range d-wave 
superconducting correlation in the interesting parameter range of J/t < 0.5. Implications of this 
result will also be discussed. 



PACS number: 74.20.-z, 71.27.+a, 74.25. Dw 

It is believed that the two-dimensional (2D) t— J model 
is a reasonable starting point to understand the phys- 
ical properties of the high temperature superconduct- 
ing (HTSC) cuprates One of the critical issues 
is whether the model has enough ingredients to quan- 
titatively explain the high T c . Here we shall report a 
numerical study to address this issue. 

The 2D t — J model on a square lattice is: 

H = -t ^ (4rC,V + h.C.) + J ^ (Si ■ Sj - -UiTlj), 
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where cj. 



C L(1 ~ n i-a)i an d denotes the near- 

est neighbors i and j. Using the experimental results of 
HTSC, we expect the physical interesting value of J/t to 
be about 0.4 and that superconductivity should occur for 
electron density n e greater than 0.7. 

The first indication about the superconductivity of the 
t — J model is to determine if the two holes would form a 
bound state by the exact diagonalization study on small 
lattices ||. However, the attractive potential among 
doped holes is not a sufficient condition for supercon- 
ductivity. The long-range pair-pair correlation should be 
a better indicator. The pair-pair correlation function is 
defined as: 



P s or d(R) = ^J2( A t A < 



l) 



(2) 



where A< = c^(c i+x i + c i-&\ ± Ci+yi ± Ci_g;). + and - 
represent s and d x 2_ y 2 pairing states respectively. N s is 
the number of sites. 

White et al. M studied the one-band Hubbard model 
and suggested that at low temperatures the pair-field sus- 
ceptibility Xd = J2nPd{R) is enhanced in the d x 2_ y 2 
channel and small for others. Other convincing results 
are from variational Monte Carlo studies pi. Although 
it did not provide quantitative values for T c , the magni- 
tude of pair-pair correlation varies with the hole density 
in a similar way as T c of HTSC. The prediction of the 



dominance of d-wave pairing, instead of s-wave, in the 
t — J model is also quite encouraging. 

However, these results are not quite consistent with a 
recent report by Zhang et al.. They studied the 2D one- 
band M and three-band Hubbard model H by using the 
constrained path Monte Carlo (CPMC) method. They 
concluded that for U/t > 4 the long-range pair-pair cor- 
relation vanishes. It becomes quite important to have a 
careful numerical study of the pairing correlation in the 
2D t — J model for larger lattices. In this paper we show 
that it is likely that there is no long-range d x 2_ y 2 pair- 
pair correlation at all for the two-dimensional t—J model 
in the physical parameter range (J/t < 0.5). 

In the variational level, the optimal state of the 2D 
t — J model for a range of parameters is the d x 2_ y 2 RVB 
trial wave function: 



RVB)=P d ]J( 



Uk 



d) I °> 



(3) 



, Afc = A(cosk x - 
/i. A is the d-wave 



with Vk/u k = A k /(e k + y /£ l + 
cosky) and e k = 2(cosk x + cosk y ) — 
superconducting order parameter and fi is the chemical 
potential. The operator Pd enforces the constraint of no 
double occupancy. We take t = 1 in this paper. This 
wave function with the form of a projected BCS wave 
function is known to be superconducting |p[. 

It is well known that the variational calculation usu- 
ally overestimates the effect of superconductivity of the 
true ground state. For the t — J model the energy is 
dominated by the nearest neighbor interaction. Hence it 
naturally leads to large A in the variational calculation. 
Heeb and Rice suggested that to examine the true 
pairing correlation, it might not be a good idea to use 
the lowest variational energy as a criterion in selecting 
the parameters of the trial function. They modified the 
function Eq.(^) with parameters that they believe can 
separate the short- and long-range parts of the correla- 
tions. They found the critical J c ~ 0.44 for the onset of 
superconducting long-range order for n e — 42/50. Here 
we modify their idea with a more systematic approach 
and provide with more rigorous analysis. 
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One of the ways to eliminate the bias induced by the 
trial wave function is to project the trial function onto 
the true ground state of the system. Recently we have 
shown that the ground state energy and many other prop- 
erties might be obtained accurately by using a particu- 
lar ground state projection method, the power-Lanczos 
method , which is a hybrid of the power method and 
the variational Lanczos method. In the power method 
it can be easily shown that if a trial wave function |^) 
is not orthogonal to the ground state, (W — H) m | \&) 
is proportional to the ground state wave function as the 
power m approaches infinity. W is an appropriately cho- 
sen constant to make the ground state energy the largest 
eigenvalue of the W — H matrix. However in practice, 
due to the Fermion sign problem the power cannot be 
too large. A better way is to improve the trial function 
| ^) by using the first order Lanczos method, i.e., to use 
| PLl) = (1 + CiH) | \fr) and then apply the power of 
W — H. C\ is a new variational parameter. The results 
described below are either calculated by | 'J), denoted by 
PLO-V, or by (W - H) m | PLl) as PLl-Pm. m = is 
the variational result, PL1-V, of the | PLl) state. 
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FIG. 1. Power dependence of the long-range pairing aver- 
age P% vg for n e = 42/50, A = 0.64, and J = 0.7. 

Clearly, if the trial function is a very good represen- 
tation of the ground state, the various correlation func- 
tions calculated in PLO-V should be almost the same as 
in PLl-Pm. On the other hand, if the trial function is 
biased in the wrong way, results of PLl-Pm will be ver 
different from PLO-V and it will correct the bias [ 
This is demonstrated in Fig.[j]. 

The averaged value, P d V9 , of the long-range part (R > 
3) of Pd{R) is plotted as a function of powers in Fig.[j] 
for A = 0.64 which gives the lowest variational energy. 
The electronic density is n e — 42/50 and J = 0.7. The 
value of P°[ vg is suppressed substantially from the VMC 
or PLO-V result when the power is increased. Clearly 
the optimized trial function has grossly overestimated the 
strength of the long range pair-pair correlation. As noted 



above, this is due to the choice of a large A to optimize 
the short-range pairing. 

Although the result in Fig.[j] is not yet converged, we 
could see that the long-range pair correlation changes 
monotonically as the wave function approaches the 
ground state. A more clear demonstration of this mono- 
tonic behavior is shown in Fig.^. 
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FIG. 2. Power dependence of the long-range pairing aver- 



;e P" V3 for several trial wave functions with different A for 
0.7. 



n e = 42/50 and J 



pavg 



is plotted as a function of powers in FigJ| for 
three different values of A: open circles are for A = 0.18, 
full circles for 0.22 and triangles for 0.26. P d V9 remains 
almost unchanged for A = 0.22. For A larger than 0.22, 
the pair correlation always decreases with power. And 
the opposite is true for A smaller than 0.22. Since P d V9 
for A = 0.22 hardly changes as the state gets closer and 
closer to the ground state, it is natural to assume that 
this is the ground state result. Moreover, if this is a 
good criterion to determine the ground state value, then 
we really only need to calculate PLO-V, PL1-V and PL1- 
P2. There is no need to go to larger powers and the 
minus sign problem is avoided. The same result would 
be obtained if we examine Pd(R) for the largest R instead 
of using P d . 

In addition to A, /z is also a variational parameter in 
the RVB wave function. Just like A which is not the 
real superconducting order parameter, /i is not the real 
chemical potential as in the simple mean field equations. 
In the discussion in the previous paragraph, fi is cho- 
sen to be consistent with the Fermi surface of the ideal 
fermi gas. For example, fi — —0.4 for n e = 42/50. If 
P d vg is the true ground state value, then just like A no 
matter what initial values of \i we use for the trial wave 
function, the converged result would remain the same. 
In other words, different sets of (A, fi) will converge to 
the same final P d V9 ■ This important consistency check 
has been verified. For example, for n e = 42/50, we ob- 
tain P e d lvg = 0.0245(4) for A J= i = 0.4 and fj, = -0.4; 
P a d V9 = 0.0238(4) for A J=1 = 0.34 and fx = -0.6; and 
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P™ 9 = 0.0238(3) for A J=1 = 0.24 and fx = -0.8. 

So far by using the ground state projection method, 
we have shown that the long-range pair correlation P^ vg 
approaches the ground-state value monotonically. This 
basic assumption is used to choose parameters to best 
represent the ground-state value of P^" 9 instead of de- 
termining variational parameters according to the lowest 
energy criterion. Unfortunately there are no exact calcu- 
lations for the 2D t — J model with large lattices to test 
our assumption. However, the method could be tested in 
the study of the pairing correlation for the 2D attractive 
and repulsive Hubbard models. 



• VMC 

V Power- 
□ Exact 



the range of values for different degenerate states of the 
IFG. The nonzero P^" 9 is obviously a finite size effect. 
Since P^" g for J = 0.5 are smaller than that of the IFG 
we can easily conclude that there is no long range d-wave 
pairing correlation for J = 0.5. 
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FIG. 4. P* vg for J ft = 0.5 and ideal fermi gas for 82 and 
122 sites. The PL0 and PL1 results of energy-optimizing A op t 
for 82 sites are shown in the insets. 



FIG. 3. Pairing correlation as a function of U for 
n e — 4/64. Full circles are results evaluated from the en- 
ergy-optimizing trial wave functions, open squares for exact 
results and open triangles are results based on our method. 

It is known that the 2D attractive Hubbard model has 
long-range s-wave pairing correlation. The on-site pair- 
ing correlation, Aj = Cj-fCjj,, for several U is shown in 
Fig.||. We consider the electron density at 4/64 which 
is solved exactly. The figure shows that the energy- 
optimizing trial wave functions (full circles) always over- 
estimate the pairing correlation in comparison with the 
exact results (open squares) . And the results obtained by 
our power-Lanczos method (open triangles) are in much 
better agreement. 

We have examined the long-range d-wave pairing cor- 
relation for the 2D repulsive Hubbard model. Our result 
agrees with Zhang et al || that the long-range pairing 
correlation is vanishingly small. 

The success for the attractive and repulsive Hubbard 
models further support our method. The method allows 
us to calculate the ground state P^ vg for large lattices. 
In Fig.[|, Pj Vg is plotted a function of electron density for 
82 and 122 sites with J = 0.5. The PL0-V and PL1-V 
results for the trial wave functions with A op t optimizing 
the variational energies are also shown in the insets of 
Fig.^j for 82 sites. It is clear that P^ v 9 pt is substantially 
reduced by PL1. And the variational values are an order 
of magnitude greater than the ground state values. 

Further, the P^ vg for J = 0.5 can be compared with 
that of the ideal fermi gas (IFG). The error bars denote 
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FIG. 5. P;" 9 for 82 sites for different density and J. 

In Fig|5| we show P^ vg for different densities and J for 
82 sites. The dotted line is the result of the IFG. It can 
be seen that P% V9 is larger than the IFG values only for 
J > 0.6, which is considerably larger than the physical 
value of J = 0.3 or 0.4. The large values observed for 
J > 0.6 could be due to the phase separation ||. For 
J < 0.5 not only P^" 9 seems to be too small to give 
large T c for the HTSC materials, the maximum P™ 9 
are at hole density Xh = 1 — x e ~ 10%. This disagrees 
with experiments which have optimal doping at ~ 
15% ~ 20%. For J > 0.6, the values of Xh are very close 
to the critical hole densities of phase separation 0], It 
is possible that for large J the P°j" 9 measured actually 
indicates electron clustering near phase separation rather 



3 



than superconductivity. 

Our result that there is no long range d-wave pairing 
correlation for J < 0.5 is actually consistent with the 
exact numerical results for the two-hole binding energy. 
It is known that binding two holes is a necessary con- 
dition for the occurrence of superconductivity. In Fig.^ 
the binding energy of two holes for various J is plotted 
as a function of the inverse of N s which is the cluster 
size. The results of 32 sites are obtained by Leung |ll]] 
and 26 sites by Poilblanc It shows that two holes 

do not bind together unless J > 0.8. Due to the different 
cluster shapes, the data are not quite on a straight line. 
However, even taking into account the deviation, the re- 
sult is still consistent with the absence of hole binding 
for J < 0.5. 
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FIG. 6. Two-hole binding energies as a function of the in- 
verse of cluster size for different J. 



bution of next-nearest neighbor hopping and next-next- 
nearest neighbor hopping, t' and t" respectively, to the 
density of states. By choosing appropriate t' and t", the 
density of states may have van Hove singularities at hole 
density around 15%. This might greatly enhance super- 
conductivity. Contrary to this belief we found that the 
superconductivity is not enhanced by adding t' term even 
though we have tuned the parameters to have the Fermi 
surface passing through the van Hove singularity at 15% 
hole density. Details will be reported elsewhere |l3| . 

Another possibility is that due to other interactions 
such as electron-phonon the effective J/ 1 might be larger 
than 0.6. However in this case, it is more likely we will 
have phase separation. The doping dependence of T c is 
also inconsistent with experiments. The intcrlayer tun- 
neling model proposed by Anderson et al certainly 
expects the absence of pairing in our pure 2D model. 
More exotic possibility might be that the true ground 
state symmetry is not d-wave but a time-reversal-brokcn 
order parameter fl5|J]~^ ] . 
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In summary, based on a simple observation that the 
long range d-wave pairing correlation changes monoton- 
ically when the state approaches the ground state, we 
assume that the ground state value can be determined 
by choosing parameters that do not optimize the varia- 
tional energy but keeps correlation value unchanged when 
we project the wave function onto the ground state. As 
a consistency check we show that different sets of val- 
ues of parameters, A and /x, produce same correlation 
value. This method has also been successfully tested for 
the 2D attractive and repulsive Hubbard models. The 
result that the long range d-wave pairing correlation is 
nonzero only for J > 0.6 is consistent with the finite size 
analysis of exact results for small clusters. It is also con- 
sistent with the result obtained in a completely different 
analysis that the phase separation boundary occurs only 
for J > 0.6. 

Since there are a number of evidences that the 2D t- 
J model is a fairly good model for HTSC, the negative 
result with respect to the d-wave pairing correlation re- 
ported above is quite disturbing. There are several pos- 
sibilities to explain this discrepancy. 

The first thing one can point to is the possible contri- 



[1 
[2 

[3] 
[4] 
[5] 
[6] 

[r] 
[s 

[9 
[10 

[11 

[12 
[13 
[14 

[is; 

[16 



P. W. Anderson, Science 235, 1196 (1987). 

F. C. Zhang, and T. M. Rice, Phys. Rev. B37, 3759 

(1988). 

D. Poilblanc et al., Phys. Rev. B 49, 12318 (1994). 
S. R. White et al., Phys. Rev. B 39, 839 (1989). 
C. Gros, Phys. Rev. B38, 931 (1988). 
Shiwei Zhang et al., Phys. Rev. Lett. 7 8, 4486 (1997). 



Shiwei Zhang et al., bond-mat/9711125| (1997). 
Y. C. Chen and T. K. Lee, Phys. Rev. B51, 6723 (1995). 

C. T. Shih et al., Phys. Rev. B57, 627 (1998). 

E. S. Heeb and T. M. Rice, Europhys. Lett 27, 673 
(1994). 

P.W. Leung, private communication. 

D. Poilblanc, Phys. Rev. B48, 3368 (1993). 

C.T. Shih, Y. C. Chen, and T. K. Lee, not published. 

S. Chakravarty et al., Science 261, 337 (1993); P. W. 

Anderson, Science 268, 1154 (1995). 

K. Krishana et al. Science 277 83 (1997). 

R. Laughlin, |oond-mat/9709004j (1997). 



4 



